Large-scale genome-wide SNP analysis reveals the rugged (and ragged) landscape of global ancestry, phylogeny, and demographic history in chicken breeds. / å¤§è§„æ¨¡å ¨åŸºå› ç»„SNPåˆ†æžæ­ç¤ºäº†é¸¡å“ç§çš„å ¨çƒç¥–å ˆã€ç§ç¾¤å‘å±•å’Œç§ç¾¤åŽ†å²çš„å¤æ‚(å’Œå¤šæ ·)çš„é—ä¼ å›¾è°±.

The worldwide chicken gene pool encompasses a remarkable, but shrinking, number of divergently selected breeds of diverse origin. This study was a large-scale genome-wide analysis of the landscape of the complex molecular architecture, genetic variability, and detailed structure among 49 populations. These populations represent a significant sample of the world's chicken breeds from Europe (Russia, Czech Republic, France, Spain, UK, etc.), Asia (China), North America (USA), and Oceania (Australia). Based on the results of breed genotyping using the Illumina 60K single nucleotide polymorphism (SNP) chip, a bioinformatic analysis was carried out. This included the calculation of heterozygosity/homozygosity statistics, inbreeding coefficients, and effective population size. It also included assessment of linkage disequilibrium and construction of phylogenetic trees. Using multidimensional scaling, principal component analysis, and ADMIXTURE-assisted global ancestry analysis, we explored the genetic structure of populations and subpopulations in each breed. An overall 49-population phylogeny analysis was also performed, and a refined evolutionary model of chicken breed formation was proposed, which included egg, meat, dual-purpose types, and ambiguous breeds. Such a large-scale survey of genetic resources in poultry farming using modern genomic methods is of great interest both from the viewpoint of a general understanding of the genetics of the domestic chicken and for the further development of genomic technologies and approaches in poultry breeding. In general, whole genome SNP genotyping of promising chicken breeds from the worldwide gene pool will promote the further development of modern genomic science as applied to poultry.

Genome-Wide Association Study Revealed Putative SNPs and Candidate Genes Associated with Growth and Meat Traits in Japanese Quail.

The search for SNPs and candidate genes that determine the manifestation of major selected traits is one crucial objective for genomic selection aimed at increasing poultry production efficiency. Here, we report a genome-wide association study (GWAS) for traits characterizing meat performance in the domestic quail. A total of 146 males from an F2 reference population resulting from crossing a fast (Japanese) and a slow (Texas White) growing breed were examined. Using the genotyping-by-sequencing technique, genomic data were obtained for 115,743 SNPs (92,618 SNPs after quality control) that were employed in this GWAS. The results identified significant SNPs associated with the following traits at 8 weeks of age: body weight (nine SNPs), daily body weight gain (eight SNPs), dressed weight (33 SNPs), and weights of breast (18 SNPs), thigh (eight SNPs), and drumstick (three SNPs). Also, 12 SNPs and five candidate genes (GNAL, DNAJC6, LEPR, SPAG9, and SLC27A4) shared associations with three or more traits. These findings are consistent with the understanding of the genetic complexity of body weight-related traits in quail. The identified SNPs and genes can be used in effective quail breeding as molecular genetic markers for growth and meat characteristics for the purpose of genetic improvement.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves.

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

Risk of Sperm Disorders and Impaired Fertility in Frozen-Thawed Bull Semen: A Genome-Wide Association Study.

Cryopreservation is a widely used method of semen conservation in animal breeding programs. This process, however, can have a detrimental effect on sperm quality, especially in terms of its morphology. The resultant sperm disorders raise the risk of reduced sperm fertilizing ability, which poses a serious threat to the long-term efficacy of livestock reproduction and breeding. Understanding the genetic factors underlying these effects is critical for maintaining sperm quality during cryopreservation, and for animal fertility in general. In this regard, we performed a genome-wide association study to identify genomic regions associated with various cryopreservation sperm abnormalities in Holstein cattle, using single nucleotide polymorphism (SNP) markers via a high-density genotyping assay. Our analysis revealed a significant association of specific SNPs and candidate genes with absence of acrosomes, damaged cell necks and tails, as well as wrinkled acrosomes and decreased motility of cryopreserved sperm. As a result, we identified candidate genes such as POU6F2, LPCAT4, DPYD, SLC39A12 and CACNB2, as well as microRNAs (bta-mir-137 and bta-mir-2420) that may play a critical role in sperm morphology and disorders. These findings provide crucial information on the molecular mechanisms underlying acrosome integrity, motility, head abnormalities and damaged cell necks and tails of sperm after cryopreservation. Further studies with larger sample sizes, genome-wide coverage and functional validation are needed to explore causal variants in more detail, thereby elucidating the mechanisms mediating these effects. Overall, our results contribute to the understanding of genetic architecture in cryopreserved semen quality and disorders in bulls, laying the foundation for improved animal reproduction and breeding.

Metabolic Rate and Egg Production in Japanese Quails Can Be Predicted by Assessing Growth Parameters of Laying Hens.

The aim of the current study was to assess the female metabolic rate and test the hypothesis that there is a relationship between the egg productivity of Japanese quails from eight breeds and their morphometric, or growth, parameters. Parameters measured were body weight (B), volume (V), and surface area (S), as well as the metabolism level expressed by the ratio S/V. The collected egg performance traits were as follows: the number of eggs produced (N), the average egg weight (W), and the total egg mass (M) (i.e., N multiplied by W). To measure the S and V values, a novel technique was developed that takes into account the similarity of the quail's body to an ellipsoid. An analysis of the relationships between productivity indicators allowed us to introduce a new index called the metabolic index, B·S/V, based on all three main growth parameters in quails. Using the values of this index, we were then able to judge indirectly the level of quails' egg productivity. We went on to assess the N, W, and M values, not only depending on the size of the bird's growth parameters but also according to the degree of their changes during quail growth. These changes were expressed as the slope angles of trend lines describing the growth process data. This approach produced more accurate results for predicting the egg productivity in terms of W and M.

Investigating adverse effects of chronic dietary exposure to herbicide glyphosate on zootechnical characteristics and clinical, biochemical and immunological blood parameters in broiler chickens.

Although the herbicide glyphosate is widely used globally and considered safe, more evidence of its adverse effects on animals and humans is accumulating. The present investigation was aimed at evaluating the impact of different glyphosate concentrations on zootechnical characteristics and clinical, biochemical and immunological blood parameters in Ross 308 broiler chickens. Four groups were employed, including untreated control and three experimental groups fed diets enriched with glyphosate at doses of 10, 20 and 100 ppm that conformed to 0.5, 1 and 5 maximum residue limits, respectively. The results showed that glyphosate is a stress factor triggering a multifaceted effect on important blood parameters (e.g., white blood cell and phagocytic counts), which was shown for the first time in the experiments involving productive meat-type poultry. It was first revealed that glyphosate-induced changes in blood parameters may be related to a negative impact on the zootechnical characteristics including the digestive tract organ development and body weight gain. The study findings suggested that exposure to glyphosate in the feedstuffs can adversely affect the physiological condition and productivity of broilers.

Genotyping-by-Sequencing Strategy for Integrating Genomic Structure, Diversity and Performance of Various Japanese Quail (Coturnix japonica) Breeds.

Traces of long-term artificial selection can be detected in genomes of domesticated birds via whole-genome screening using single-nucleotide polymorphism (SNP) markers. This study thus examined putative genomic regions under selection that are relevant to the development history, divergence and phylogeny among Japanese quails of various breeds and utility types. We sampled 99 birds from eight breeds (11% of the global gene pool) of egg (Japanese, English White, English Black, Tuxedo and Manchurian Golden), meat (Texas White and Pharaoh) and dual-purpose (Estonian) types. The genotyping-by-sequencing analysis was performed for the first time in domestic quails, providing 62,935 SNPs. Using principal component analysis, Neighbor-Net and Admixture algorithms, the studied breeds were characterized according to their genomic architecture, ancestry and direction of selective breeding. Japanese and Pharaoh breeds had the smallest number and length of homozygous segments indicating a lower selective pressure. Tuxedo and Texas White breeds showed the highest values of these indicators and genomic inbreeding suggesting a greater homozygosity. We revealed evidence for the integration of genomic and performance data, and our findings are applicable for elucidating the history of creation and genomic variability in quail breeds that, in turn, will be useful for future breeding improvement strategies.

Probiotics as an alternative to antibiotics in modulating the intestinal microbiota and performance of broiler chickens.

AIMS: Gut bacteria play an important role in poultry nutrition and the immune defense system. Changes in the intestinal microbiome affect the physiological state, metabolism, and innate immunity of poultry. The present study aimed to characterize age-related changes in the gastrointestinal tract microflora in broiler chickens, depending on supplementation of the diet with the in-feed antibiotic Stafac® 110 and a Bacillus subtilis strain-based probiotic. METHODS AND RESULTS: In this regard, a comprehensive analysis of the taxonomic structure of the microbial community in the gastrointestinal tract (GIT) of broiler chickens was carried out using a molecular genetic technique of the terminal-restriction fragment length polymorphism analysis and taking into account age dynamics and feeding treatment. A beneficial effect on the microbiological composition and body weight of broilers was observed when using the antibiotic and probiotic in compound feeds. Different bacterial communities were revealed in the duodenum and cecum, and their positive impact on broiler growth was established. The results obtained shed light on the formation of GIT microflora of broiler chickens during the growing period and its changes in response to the use of the antibiotic and the probiotic. CONCLUSIONS: We suggest that the implementation of the tested in-feed antibiotic and probiotic can be beneficial in regulating the intestinal microflora microbiological processes in the GIT and improving the feeding efficiency and productivity of broiler chickens.

A novel model for eggs like pears: How to quantify them geometrically with a mathematical formula of two parameters?

Of the variety of bird egg shapes, perhaps the most fascinating and unusual are pyriform (pear-shaped, or conical) eggs. Among oologists, there is still no consensus on what exactly caused this evolutionary and ecological adaptation. To address this, our research was aimed to develop a minimalistic mathematical model for an accurate description of the pyriform egg contour, using the minimum number of measurements. As such, egg length (L) and its maximum breadth (B) were found to be an optimal set of parameters that were easy enough to measure with the required accuracy. We tested four analytical premises that can be used for successful pyriform egg shape modelling. To validate these four model premises, images of pyriform eggs characteristic of 32 species were used. As a result, we derived a novel mathematical dependence that we called the 'pyriform model with two parameters'. Based on this model, it is feasible to geometrically reconstruct any pyriform egg profile under study using only two egg measurements, i.e., L and B. Since pyriform eggs are characteristic of wild bird species that are most often investigated in the field, the measurement of only two parameters minimizes the time spent and, accordingly, the stress factor on the animals. The least error estimate for the new model was 3.9%, which turned out to be even more accurate than that of the previously developed model with three parameters.

A novel approach to egg and math: Improved geometrical standardization of any avian egg profile.

Developing a geometric formulation of any biological object has a number of justifications and applications. Recently, we developed a universal geometric figure for describing a bird's egg in any of the possible basic shapes: spherical, ellipsoidal, ovoid, and pyriform. The formulation proved widely applicable but had a number of drawbacks, including a very obvious join between two parts of the egg. To correct this, we developed the Main Axiom of the universal mathematical formula. This essentially involved making the ordinate of the extremum of the function correspond to half the maximum egg breadth (B), and the abscissa to the reciprocal of the parameter w that reflects the shift of the vertical axis to its coincidence with B. This, in turn, helped us develop a new, simplified mathematical model without a nonbiological join. Experimental verification was performed to confirm the adequacy of the new geometric figure. It accurately described actual avian eggs of various shapes more closely than our previous model. To the best of our knowledge, our new, simplified equation can be applied as a standard for any bird egg that exists in nature. As a rather simple equation, it can be used in a broad range of applications.

Whole Genome Screening Procures a Holistic Hold of the Russian Chicken Gene Pool Heritage and Demographic History.

A study for genomic variation that may reflect putative selective signaling and be associated with economically important traits is instrumental for obtaining information about demographic and selection history in domestic animal species and populations. A rich variety of the Russian chicken gene pool breeds warrants a further detailed study. Specifically, their genomic features can derive implications from their genome architecture and selective footprints for their subsequent breeding and practical efficient exploitation. In the present work, whole genome genotyping of 19 chicken breeds (20 populations with up to 71 samples each) was performed using the Chicken 50 K BeadChip DNA chip. The studied breed sample included six native Russian breeds of chickens developed in the 17th-19th centuries, as well as eight Russian chicken breeds, including the Russian White (RW), created in the 20th century on the basis of improving local chickens using breeds of foreign selection. Five specialized foreign breeds of chickens, including the White Leghorn (WL), were used along with other breeds representing the Russian gene pool. The characteristics of the genetic diversity and phylogenetic relationships of the native breeds of chickens were represented in comparison with foreign breeds. It was established that the studied native breeds demonstrate their own genetic structure that distinguishes them from foreign breeds, and from each other. For example, we previously made an assumption on what could cause the differences between two RW populations, RW1 and RW2. From the data obtained here, it was verified that WL was additionally crossed to RW2, unlike RW1. Thus, inherently, RW1 is a purer population of this improved Russian breed. A significant contribution of the gene pool of native breeds to the global genetic diversity of chickens was shown. In general, based on the results of a multilateral survey of this sample of breeds, it can be concluded that phylogenetic relationships based on their genetic structure and variability robustly reflect the known, previously postulated and newly discovered patterns of evolution of native chickens. The results herein presented will aid selection and breeding work using this gene pool.

Integrative comparative analysis of avian chromosome evolution by in-silico mapping of the gene ontology of homologous synteny blocks and evolutionary breakpoint regions.

Avian chromosomes undergo more intra- than interchromosomal rearrangements, which either induce or are associated with genome variations among birds. Evolving from a common ancestor with a karyotype not dissimilar from modern chicken, two evolutionary elements characterize evolutionary change: homologous synteny blocks (HSBs) constitute common conserved parts at the sequence level, while evolutionary breakpoint regions (EBRs) occur between HSBs, defining the points where rearrangement occurred. Understanding the link between the structural organization and functionality of HSBs and EBRs provides insight into the mechanistic basis of chromosomal change. Previously, we identified gene ontology (GO) terms associated with both; however, here we revisit our analyses in light of newly developed bioinformatic algorithms and the chicken genome assembly galGal6. We aligned genomes available for six birds and one lizard species, identifying 630 HSBs and 19 EBRs. We demonstrate that HSBs hold vast functionality expressed by GO terms that have been largely conserved through evolution. Particularly, we found that genes within microchromosomal HSBs had specific functionalities relevant to neurons, RNA, cellular transport and embryonic development, and other associations. Our findings suggest that microchromosomes may have conserved throughout evolution due to the specificity of GO terms within their HSBs. The detected EBRs included those found in the genome of the anole lizard, meaning they were shared by all saurian descendants, with others being unique to avian lineages. Our estimate of gene richness in HSBs supported the fact that microchromosomes contain twice as many genes as macrochromosomes.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds.

BACKGROUND: The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes. RESULTS: Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) FST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF). CONCLUSION: Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Features of Fractal Conformity and Bioconsolidation in the Early Myogenesis Gene Expression and Their Relationship to the Genetic Diversity of Chicken Breeds.

Elements of fractal analysis are widely used in scientific research, including several biological disciplines. In this study, we hypothesized that chicken breed biodiversity manifests not only at the phenotypic level, but also at the genetic-system level in terms of different profiles of fractal conformity and bioconsolidation in the early myogenesis gene expression. To demonstrate this effect, we developed two mathematical models that describe the fractal nature of the expression of seven key genes in the embryonic breast and thigh muscles in eight breeds of meat, dual purpose, egg and game types. In the first model, we produced breed-specific coefficients of gene expression conformity in each muscle type using the slopes of regression dependencies, as well as an integral myogenesis gene expression index (MGEI). Additionally, breed fractal dimensions and integral myogenesis gene expression fractal dimension index (MGEFDI) were determined. The second gene expression model was based on plotting fractal portraits and calculating indices of fractal bioconsolidation. The bioconsolidation index of myogenesis gene expression correlated with the chick growth rate and nitric oxide (NO) oxidation rate. The proposed fractal models were instrumental in interpreting the genetic diversity of chickens at the level of gene expression for early myogenesis, NO metabolism and the postnatal growth of chicks.

Disentangling clustering configuration intricacies for divergently selected chicken breeds.

Divergently selected chicken breeds are of great interest not only from an economic point of view, but also in terms of sustaining diversity of the global poultry gene pool. In this regard, it is essential to evaluate the classification (clustering) of varied chicken breeds using methods and models based on phenotypic and genotypic breed differences. It is also important to implement new mathematical indicators and approaches. Accordingly, we set the objectives to test and improve clustering algorithms and models to discriminate between various chicken breeds. A representative portion of the global chicken gene pool including 39 different breeds was examined in terms of an integral performance index, i.e., specific egg mass yield relative to body weight of females. The generated dataset was evaluated within the traditional, phenotypic and genotypic classification/clustering models using the k-means method, inflection points clustering, and admixture analysis. The latter embraced SNP genotype datasets including a specific one focused on the performance-associated NCAPG-LCORL locus. The k-means and inflection points analyses showed certain discrepancies between the tested models/submodels and flaws in the produced cluster configurations. On the other hand, 11 core breeds were identified that were shared between the examined models and demonstrated more adequate clustering and admixture patterns. These findings will lay the foundation for future research to improve methods for clustering as well as genome- and phenome-wide association/mediation analyses.

Landmark native breed of the Orenburg goats: progress in its breeding and genetics and future prospects.

This paper reviews information about a unique and iconic breed of the Orenburg Oblast, the homeland and the only place where the best herds of Orenburg down-hair goats in Russia are concentrated. Three types of these small ruminant animals are widespread on the territory of the region: Orenburg purebred gray goats, Orenburg purebred white goats, as well as crossbred white goats of F1 White Don × White Orenburg. Currently, at the farms of the Orenburg region, animals are selected according to their phenotype, with selected traits being color, weight and length of down hair. In recent years, the Orenburg goat breed has become an object of genetic research using various marker systems including immunogenetic, microsatellite, mtDNA and SNP markers. Overall, these studies evidence about the uniqueness of the allele pool in the landmark native breed of the Orenburg goats, which is a complex dynamic genetic system, prioritizing its further in-depth genome research and breeding applications.

Towards a Genetic Linkage Map of the California Condor, an Endangered New World Vulture Species.

The development of a linkage map is an important component for promoting genetic and genomic studies in California condors, an endangered New World vulture species. Using a set of designed anonymous microsatellite markers, we genotyped a reference condor population involving 121 individuals. After marker validation and genotype filtering, the genetic linkage analysis was performed using 123 microsatellite loci. This resulted in the identification of 15 linkage groups/subgroups that formed a first-generation condor genetic map, while no markers linked to a lethal chondrodystrophy mutation were found. A panel of polymorphic markers that is instrumental in molecular parentage diagnostics and other genetic studies in the California condor was selected. Further condor conservation genomics research will be focused on updating the linkage map and integrating it with cytogenetic and BAC-based physical maps and ultimately with the genome sequence assembly.

Egg-inspired engineering in the design of thin-walled shelled vessels: a theoretical approach for shell strength.

A novel subdiscipline of bionics is emerging in the form of 'egg-inspired engineering' through the use of egg-shaped ovoids as thin-walled tanks and building structures. Hügelschäffer's and Narushin's models of egg geometry are highly applicable within this proposed subdiscipline. Here we conducted a comparative analysis between the two models with respect to some of the most important egg parameters. These included contents volume, shell volume, and the location of the neutral axis along the shell thickness. As a first step, theoretical studies using the Narushin's model were carried out due to the lack (or limited amount) of data on the geometric relationships of parameters and available calculation formulae. Considering experimental data accumulated in the engineering and construction industries, we postulate a hypothesis that there is a correlation between location of the neutral axis and the strength of the walls in the egg-shaped structure. We suggest that the use of Narushin's model is preferable to Hügelschäffer's model for designing thin-walled shelled vessels and egg-shaped building structures. This is due to its relative simplicity (because of the requirement for only two initial parameters in the basic equation), optimal geometry in terms of material costs per unit of internal capacity, and effective prerequisites for shell strength characteristics.

Comparative Analysis of Molecular RFLP and SNP Markers in Assessing and Understanding the Genetic Diversity of Various Chicken Breeds.

Monitoring the genetic diversity of small populations is important with respect to conserving rare and valuable chicken breeds, as well as discovery and innovation in germplasm research and application. Restriction fragment length polymorphisms (RFLPs), the molecular markers that underlie multilocus DNA fingerprinting (MLDF), have historically been employed for this purpose, but over the past two decades, there has been an irreversible shift toward high-throughput single-nucleotide polymorphisms (SNPs). In this study, we conducted a comparative analysis of archived MLDF results and new data from whole-genome SNP genotyping (SNPg) among 18 divergently selected breeds representing a large sample of the world gene pool. As a result, we obtained data that fit the general concept of the phylogenetic distribution of the studied breeds and compared them with RFLP and SNP markers. RFLPs were found to be useful markers for retrospective assessment of changes in the genetic architecture and variability underlying the phenotypic variation in chicken populations, especially when samples from previous generations used for MLDF are unavailable for SNPg. These results can facilitate further research necessary to assess the possibility of extrapolating previous MLDF results to study the long-term dynamics of genetic diversity in various small chicken germplasm populations over time. In general, the whole-genome characterization of populations and breeds by multiple SNP loci will further form the basis for the development and implementation of genomic selection with the aim of effective use of the genetic potential of the domestic gene pool in the poultry industry.

Measurement of the neutral axis in avian eggshells reveals which species conform to the golden ratio.

Avian eggs represent a striking evolutionary adaptation for which shell thickness is crucial. An understudied eggshell property includes the neutral axis, a line that is drawn through any bent structure and whose precise location is characterized by the k-factor. Previous studies have established that, for chicken eggs, mean k corresponds to the golden ratio (Φ = 1.618, or 0.618 in its reciprocal form). We hypothesized whether such an arrangement of the neutral axis conforms to the eggshell of any bird or only to eggshells with a certain set of geometric parameters. Implementing a suite of innovative methodological approaches, we investigated variations in k of 435 avian species, exploring which correspond to Φ. We found that mean k is highly variable among birds and does not always conform to Φ, being much lower in spherical and ellipsoid eggs and higher in pyriform eggs. While 21 species had k values within 0.618 ± 0.02 (including four falcon species) and the Falconinae subfamily (six species) revealed a mean of 0.618, it is predominantly domesticated species (chicken, ducks, and geese) that lay eggs whose neutral axis corresponds to the golden ratio. Thus, the study of the mathematical secrets of the eggshell related to the golden ratio of its neutral axis suggests its species-specific signatures in birds.